Spin parity and broken symmetry in finite spin-1/2 chains with frustrated exchange: 

quantum transition from high to low spin 
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Exact diagonalization of finite spin-1/2 ciiains witli periodic boundary conditions is applied to 
the ground state (gs) of chains with ferromagnetic (F) exchange Ji < between first neighbors, 
antiferromagnetic (AF) exchange J2 — aj\ > between second neighbors, and axial anisotropy 
< A < 1. In zero field, the gs is in the Sz = Q sector for the relevant parameters and is doubly 
degenerate at multiple points 7^, = (Qm, Am) in the a, A plane. Degeneracy under inversion at sites 
or spin parity or both leads, respectively, to a bond order wave (BOW), to staggered magnetization 
or to vector chiral (VC) order. Exact results up to A'^ = 28 spins directly yield order parameters 
and spin correlation functions whose weak N dependencies allow inferences about infinite chains. 
The high-spin gs at J2 = changes discontinuously at 71 = ( — 1/4, 1) to a singlet in the isotropic 
(A = 1) chain. The transition from high to low spin S{a, A) is continuous for A < As = 0.95 ±0.01 
on the degeneracy line q:i(A). The gs has staggered magnetization between Aa = 0.72 and As, 
and a BOW for A < Aa. When both inversion and spin parity are reversed at 7™, the correlation 
functions C(p) for spins separated by p sites are identical. C(p) minima are shifted by 7r/2 from 
the minima of VC order parameters at separation p, consistent with right and left-handed helices 
along the z axis and spins in the xy plane. Degenerate gs of finite chains are related to quantum 
phase diagrams of extended a, A chains, with good agreement for order parameters along the line 
ai(A). Degenerate gs limit a VC phase to intermediate a and A where ^(a. A) varies rapidly but 
continuously, in contrast to many-body treatments in which VC phases extend over a larger range 
of parameters. 
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I. INTRODUCTION 

One-dimensional (ID) spin chains have been exten- 
sively studied over the years both experimentally and 
theoretically. Spin chains are good approximations 
to the magnetism of diverse inorganic and organic 
crystals. They are simple many-body quantum systems, 
well suited for computational studies with some 
exactly known properties and rich ground-state (gs) 
phase diagrams. Solid-state studies focus broadly on 
magnetic properties, instabilities and phase transitions 
that limit ID behavior at low temperature. Theoretical 
interest extends to quantum phases with different gs in 
parameter space. In this paper we consider the gs prop- 
erties of spin-1/2 chains, Eq. [T] below, with isotropic 
or axially anisotropic exchange Ji between nearest 
neighbors and antiferromagnetic (AF) exchange J2 > 
between second neighbors that ensures frustration for 
either sign of Ji. The spin-Peierls system CuGeOa 
illustrates spin-1/2 chains of Cu(II) ions with Ji > 0. 
The isotropic AF/AF chain is a spin liquid up to [3] 
a = J2/J1 < 0.2411 and its exact gs is a bond order 
wave (BOW) at a = 1/2, the Majumdar- Ghosh point 4]. 

Spin-1/2 Cu(II) chains with ferromagnetic (F) Ji < 
have recently been identified in cupric oxides |5l4l3| with 
estimated a's ranging frorn [l3ll a w in Ba3Cu3ln40i2 
or Ba3Cu3Sc40i2 to [H, [g a w -0.5 in LiCuSb04, 



LiCuZr04, LiCuV04 and LiCu202. The isotropic F/AF 
chain has a ferromagnetic gs for a — J2I J\ 1^ etc — —1/4, 
a singlet gs for a < ac, and exact degeneracy at the 
quantum critical point Oc as shown by Hamada et al.fisj 
Vector chiral (VC) and multipolar phases have been of 
special interest [ij, [16142 Oj . Hikihara et al. [1^ discuss 
the phase diagram of the isotropic F/AF chain in a static 
magnetic field. Furukawa et al. [l7[ and Sirker [ij], 
among others [2l|, have studied the axially anisotropic 
F/AF chain in zero field, the model considered in this 
paper. The limit Ji = decouples the system into two 
AF chains as sketched in Fig. [T]for a zigzag chain. 

The Hamiltonian of the anisotropic F/AF chain with 
periodic boundary conditions (PBC) and spin-1/2 sites 
is 

N 

H{a,A) = Ji y^X^p ■ Sp+i + aSp ■ Sp+2) 
p=i 

+ {A-i){s;s;+, + as;s;+2) (i) 

Ji = — 1 sets the energy scale. The model has two 
parameters, the frustration ratio a — J2/J1 < and 
axial anisotropy < A < 1. Total spin S is conserved in 
the isotropic limit, A = 1, but only Sz is a good quantum 
number otherwise. The gs is always in the Sz — sector 
for the parameters of interest in this work, and Sz = 
basis states are products of A^/2 spins a and N/2 spins 
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/?. The gs transforms as P = ±1 under the spin parity 
operator P that reverses all spins and as Q — ±1 
under inversion at sites, which corresponds to reflection 
through sites p, p + N/2 in finite systems. H{a, A) also 
has Cn translational symmetry and inversion symmetry 
midway between sites. 

The gs phase diagram of H(a,A) has been studied by 
many-body methods P, [13, [l6[ - [23j such a field theory, 
perturbation theory and density matrix renormalization 
group (DMRG) that so far agree only in part. Broken 
symmetry is expected and found for ranges of a and 
A, with multiple exotic phases near the quantum 
critical point ac — —1/4, A = 1. Excitation energies 
are small, of order for N spins and exceptionally 
small according to one field theory [24]. One challenge 
is to distinguish between strictly degenerate gs that 
indicate broken symmetry and nondegenerate gs with 
tiny excitation energies, as discussed carefully by Afileck 
and Lieb in ID spin chains j25j . 

We adopt a different approach that bears directly on 
broken symmetry and order parameters but only indi- 
rectly on the phase diagram. We solve H{a, A) exactly 
for finite N npto N = 28. Since exact eigenstates respect 
all symmetries, broken symmetry requires degenerate 
gs that in turn yield order parameters. On the other 
hand, phase boundaries from finite-size calculations 
are based on excited-state crossovers, or degeneracy, as 
discussed for the BOW phase of the isotropic AF/AF 
chain 0, ^ or of extended Hubbard models [13, ^ 
Order parameters become nonzero at the boundary and 
increase in the broken-symmetry phase. As shown in 
Section II, H{a, A) has doubly degenerate gs at multiple 
parameter values = (am, A^), = 1,2,..., where 
the gs symmetry changes. 

There are three kinds of gs degeneracy. When 7™ cor- 
responds to Ci = ±1 or P = ±1, the linear combinations 



|7m) = (l7m,l)±l7m,-l))/V2 



(2) 



have broken inversion or spin parity symmetry and differ- 
ent order parameters. When both d and P are reversed 
at either (1,1) and (-1,-1) or (1,-1) and (-1,1), the 
linear combinations in Eq. [5] have CiP — ±1. There 
are finite-size gaps to all other eigenstates, sometimes re- 
markably small gaps. Broken Q symmetry leads as usual 
to a BOW or dimer phase; broken parity P is associ- 
ated with staggered magnetization or Neel order; broken 
Ci and P symmetry yields VC order. Dimer, Neel and 
VC phases are all possibilities for H{a, A) in zero field 
[3 [17, 19]. Exact diagonalization, albeit limited to dis- 
crete 7m, makes it possible to compute order parameters 
as well as spin correlations functions or other properties. 
When systematic variations with N are found, extrapola- 
tion gives accurate but not exact information about the 
infinite chain. 




FIG. 1: Schematic representation of a spin-1/2 chain with 
exchange Ji and J2 between first and second neighbors. 



Increasing J2 (increasing —a) induces a quantum tran- 
sition from a high to low spin. The isotropic chain has 
a first order transition with discontinuous S" at 71 — 
(—1/4, 1). Sz rather than S is conserved in systems with 
axial anisotropy A < 1. The normalized spin per site 
S{a, A) is the gs expectation value 



2S'(a,A) 



2{a, A\S^\a, Ay/^ 
[iV(Ar + 2)]i/2 



< 1 



(3) 



S{a, A) is shown in Section III to be discontinuous for 
small anisotropy A > A^ ~ 0.95, continuous for A < 
Ab- There is finite Neel order between A^ = 0.75 and 
Ab, finite BOW or dimer order for A < A^. The spin 
transition requires exact gs and has apparently not been 
recognized previously in anisotropic chains. Spin correla- 
tion function show spiral order at 7m when Cj symmetry 
is broken, as in isotropic AF/AF chains [2a,[23|- We find 
identical spin correlations functions at parameter values 
7m when both Ci and P are broken, and interpret these 
results as right and left-handed helices along the unique 
axis with spins in the xy plane. We compare finite-size 
results with previous theory and find considerable agree- 
ment as well as occasional disagreement. Our results 
limit VC degeneracy to small ranges of parameters a, A 
in which S{a, A) varies rapidly but continuously. More 
extended VC phases have been inferred by other methods 

[13, [a. 



II. DEGENERACY AND BROKEN SYMMETRY 

We summarize some properties of the anisotropic 
F/AF spin chain, Eq. [TJ before presenting numerical 
results for even N with PBC. The gs is ferromagnetic 
at J2 = (a = 0) with magnetization in the xy 
plane for A < 1. Increasing J2 > (increasing —a) 
induces a quantum transition to a low-spin state. The 
ferromagnetic and singlet gs of the isotropic chain 
are degenerate at Uc = —1/4 with energy per site 
eo = —3/16. The exact result for the extended system is 
the first degeneracy 71 = (—1/4,1) for finite N, where 
S changes from N/2 to 0. Any singlet with Sz ~ can 
be represented [s^] as linear combinations of N/2 paired 
spins {ai/3j — /3iaj)/^ whose phase is fixed by choosing 
site i < j. The parity is P = (—1)^/^ since reversing 
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all spins gives a phase factor of -1 for every singlet 
pair. The gs linear combination at ac is the uniformly 
distributed resonating valence bond solid [l5l |. 

The degeneracy 71 = ai{A) between high and low 
spin shifts to ai(A) < ac = —1/4. The unit step 
function 2S{a,l) at ac = —1/4 decreases rapidly for 
A < 1 and 2S{a, A) is continuous for A < 0.95. In 
the limit A = of extreme anisotropy, the crossover 
to low spin occurs at ai(0) = —1/2 for even N. The 
degeneracy 71 = (—1/2,0) is the F/AF version of the 
Majumdar- Ghosh point of the isotropic AF/AF chain [4|. 
Here the exact gs is a product of triplets with Sz = 0, 
{tti^j + aiPj)/V2, either with i — 2n — 1, j = 2n ot 
i = 2n, j = 2n + 1. There are N/2 triplets and Eq. [3] 
gives 



5(-l/2,0) = l/v/(iVT2) 



(4) 



The exact result shows that the extended system has 
S{a,A) = for —a > 1/2 over the entire range of A. 
The singlet gs at A = 1 becomes a linear combination of 
states with S'^ ^ N in anisotropic chains, and the spin 
per site goes as N^^/"^ on the low-spin side. 

Another relevant limit is Ji — 0, when H{a, A) 
decouples into chains with anisotropic J2 > between 
nearest neighbors as sketched in Fig. ^ N = An systems 
decouple into ID chains of 2n sites that correspond to 
the XXZ Heisenberg spin-1/2 antifcrromagnet. The gs 
of each chain has Sz — ^ and there are no correlations 
between spins in different chains. On the other hand, 
N — 4n + 2 systems decouple into two radicals with 
Sz = ±1/2 whose gs remain entangled even at Ji = 0. 
We find 4n, An + 2 effects in isotropic chains around 
a sa —1/2, but none at 71 — ai(A). 

We obtain the lowest energy of H(a,A), Eq. [1] for 
even N and PBC in four sectors with Sz = 0, Ci = ±1 
and P — ±1. The absolute gs is degenerate at points 
7m = {arm Am) that are found to two or three significant 
figures. Preliminary scans for small on a rougher grid 
identifies parameter ranges with symmetry crossovers. 
Fig. [2] shows the 7m in the —a, A plane for A^ = 20. 
Open circles indicate Ci degeneracy, closed circles P 
degeneracy and stars VC degeneracy. The points ai at 
A = 1 and are exact, independent of size, and the 
line 71 — ai(A) hardly depends on A^. All degeneracy 
for A < 0.7 is associated with inversion symmetry. 
Staggered magnetization or VC order is limited to 
intermediate anisotropy in Fig. [21 The intermediate 
region is expanded in Fig. [3] for A = 24. 

Table [T] lists the n — 1 values of (—am, 1), in > 2, for 
isotropic N ^ An chains, all with Ci — ±1. Anisotropic 
A^ = 4n chains with A < 0.7 also have n degenerate 
points am with Ci = ±1. Broken-symmetry gs are given 
by Eq. [H When 7™ corresponds to = ±1, the BOW 
or dimer amplitude B(^rn) is 




0.4 0:5" 

-a=J2/|J^i 



FIG. 2: Parameter points (q. A) for a doubly degenerate 
ground state of H{a,A), Eq. [1] with A'^ = 20 spins. Open 
circles indicate inversion degeneracy d = ±1, closed circles 
indicate spin parity degeneracy P = ±1, and stars indicate 
vector chiral degeneracy dP — ±1. The line ai(A) marks 
degeneracy between high and low spin. 




< 8 - High Spin • 

A 

0.75- Y^=a^(A)i»»* *) 



0.3 0.4 
-a=J„/|J. 



FIG. 3: Same as Fig. [2]for A = 24 at intermediate anisotropy 
A and frustration a. 



S(7m) = (7rn; -l|(5l • 5*2 - ^2 ' ^3)|7m; l>/2 (5) 

The order parameter is the off-diagonal matrix element, 
±B for the two gs, and we used translational symmetry. 
Only BOWs are realized in isotropic (A — 1) systems 
whose B{am, 1) are shown in Fig. Ufor both N — An and 
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TABLE I: Degenerate ground states at (~am,l), m > 2, of 
isotropic spin chains, Eq. [U with A'' sites and A = 1 
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FIG. 4: (Color online) Discrete order parameter B{am), Eq. 
5, as a function of frustration a in isotropic (A = 1) chains 
of N sites; continuous B{a) from DMRG. 



Kz(p,7m) = (7m; -1, -IK'S*! X S'i+p)^|7„i; 1, 1) 

2 -MiS^Si+p - S'rS'i++p)|7m; 1, 1) 

(7) 

The z component of the vector product is finite for axial 
anisotropy, and Eq. [7] is for degeneracy between (1,1) 
and (-1,-1). The matrix element has (-1,1) and (1,-1) 
for the other way of reversing both d and P. We paid 
close attention to degenerate points 7™ and found only 
two- fold degeneracy. In a few cases, the gap to the first 
excited state is tiny, of the order of 10~^, far less than 
1/A^. Double degeneracy implies one broken symmetry 
in finite chains for any (a. A) in the range —a > 1/4, 
< A < 1. 

The order parameters 5(71), its transverse part B±, 
and Mst (71) are shown in Fig. [5] as a function of 
anisotropy along 71 — q;i(A). The size dependence is 
negligible, and A = is exact, with B = 1/8, B± = 1/4. 
As mentioned above, finite i?(7i) or Mst(7i) implies that 
the phase boundary where the order parameter becomes 
nonzero is on the high-spin side of 71 , but deviations from 
71 are only significant for large order parameters. The 
largest difference occurs at A = 0, where ai = —0.50 
and excited-state crossovers give a = —0.325. The 
dimer/Neel boundary in Fig. [S]is Aa{N) = 0.722 with 
aA = -0.317 for N = 18, 20 and 22, while the Neel/first 
order boundary is As (TV) = 0.941 with as = -0.2674. 
The staggered magnetization increases with anisotropy in 
the Neel phase to 0.050 at Aa- B{ai) is small (< 0.03) 



4n + 2. B{a, 1) has a broad maximum around a ~ —0.5 
and a modest size dependence. DMRG results in Fig. 
|4] are based on an algorithm [3l| that adds four spins 
per step, as needed for accuracy at —a > 0.5, and gives 
fragments with an even number of sites at each step. The 
AF/AF chain has a broad maximum 28] 5(0.4, 1) = 0.40 
that is 10-fold larger. The z and transverse parts of 
B{am) are found separately in anisotropic systems 
with A < 1. Most 7m in Figs. [5] and |3] refer to de- 
generacy under inversion, and i?(— 1/2, 0) = 1/8 is exact. 

When 7m corresponds to P — ±1, the amplitude 
^stilm) > of the staggered magnetization is 




M.t(7m) = 2(7m;-l|(^f - ^2^)17™; 1|) (6) 

Similarly, when both Ci and P are reversed at 7„i, the 
VC order parameters ±Kz{p'T^m) for spins p sites apart 
are 



FIG. 5: (Color online) Order parameters along the degener- 
acy line ai(A) between high and low spin. B and B± are 
Eq. [5] and its transverse part, and Mst is Eq. (6] The spin 
transition of the infinite chain is first order for A > As and 
has vanishing order parameters. 
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and decreases with N near in the first-order region 
where we expect i? = in the extended system. For 
comparison, the infinite time-evolving block decimation 
(iTEBD) algorithm |T3l gives = 0.72, = -0.320 
and As = 0.93, as = -0.272 in Fig. 4 of ref. 17a; the 
maximum Neel amplitude is 0.055 and 2B± (called -D^fa) 
= 0.35 at A = 0.65 where we find 2B± = 0.309 at iV = 
16 and 0.307 at TV = 22. There is remarkably close agree- 
ment between two completely different calculations. We 
disagree near A = 1 where iTEBD returns small finite B 
rather than B = 0. 



III. SPIN TRANSITION AND CORRELATION 
FUNCTIONS 

We evaluate the spin per site 2S{a,A), Eq. [3l for 
slightly anisotropic systems along the degeneracy 71 — 
ai(A) between high {S+) and low spin {S-). The size 
dependence of <S'+(A) — ^-{A) is shown in Fig. [HI up to 
= 28. The step function of isotropic chains is quickly 
lost with increasing N. Although 1/A^ behavior is ap- 
proximate at best, the gap has vanished by A = 0.94 
and the extended system has a continuous transition 
for A < Ai3 0.95 ±0.01. The S{a,A) discontinuity 
at ai(A) clearly decreases very rapidly with anisotropy. 
The precise value of A^ is less important than recog- 
nizing a first order quantum transition for A > As- Or- 
thogonality in S then ensures vanishing order parameters 
for an infinite chain with a first order transition in Fig. [S] 
As an indication of internal consistency, we note that al- 
most identical values of Ab are inferred from closing the 
5+ (A) — ^-(A) discontinuity and from the Neel/first- 
order boundary. 

We define gs correlation functions for spins p sites 
apart using the translational symmetry of H{a, A) 



C(P,7) = (71(^^1 -^1+^)171) 



(8) 



The z and transverse components are computed sepa- 
rately and are not simply related in anisotropic systems. 
For 5* = 1/2 chains, we have 



p=i 



{7\Sl\l)/N = i+^C^(p,7) 



(9) 



p=i 



Since the gs is in the Sz = sector, finite S{a, A) is 
due to transverse components. The high-spin regime 
with —7 < ai(A) has C{p) > for all p, or simply 
C{p) =:l/4atA = l. The low-spin regime necessarily 
has C{p) < for some p, and PBC implies an even 
number of sign changes as a function of p. 
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FIG. 6: (Color online) Discontinuity in S{a, A), Eq. [3l at the 
degeneracy ai(A) between high and low spin in system of N 
spins with axial anisotropy A up to N = 28. 
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FIG. 7: (Color online) Spin correlation functions C{6), Eq. |8] 
with 6 = 2np/N , of isotropic chains H{a, 1), Eq. [T] Upper 
panel: 02 up to = 28 with 2r changes of sign; the r = 1 and 
2 lines have B — 0.276, A = —0.0426. Lower panel: closed 
symbols, Q5 for A'' = 20 with 10 sign changes; open symbols, 
decoupled chains with Ji — and C{p) — for odd p. 



The upper panel of Fig. [7] shows C(p, 02) in isotropic 
chains up to = 28 as a function oi 6 = 2ttp/N 
with p — 1,2,...N — 1. At A = 1, there are two sign 
changes in the interval [—1/4, Q!2] and four sign changes 
in [a2,a3]. C(p, 02) is double valued. The systematic 
size dependence allows inferences about the extended 
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system. The line in Fig. [7^ has small deviations 
from sinusoidal due to C(0) = 3/4, whose contribution 
decreases as 1/iV. C{p,j) changes sign 2n times for 
—7 > In, as shown in Fig. [7]for — 20, lower panel, 
for isotropic chains at 75 = 0:5(1). For comparison, we 
include the limit Ji = 0, when C{p) = for odd p and 
has alternating sign for even p. The Ji = calculation is 
for a longer chain that shows decreasing AF correlation 
with increasing p. Although a„(l) w —1/2 in Table H] is 
far from the Ji — limit (a — >■ —00) of decoupled AF 
chains, the spin correlation functions are already similar. 
They change sign at most 2n times and account fully for 
the gs degeneracy of isotropic N = An chains. Isotropic 
chains, either F/AF or AF/AF, are limited to BOW or 
dimer phases, in clear disagreement with the VC phase 
in zero field in Fig. la of ref 19 or, in a smaller range, in 
Fig. 2b of ref 17b. 

C(p;7„i) of anisotropic chains are not sinusoidal but 
still have 2m — 2 and 2m sign changes when 7^ marks 
broken Ci or P symmetry. The C(p, 72) in Fig. [S] are 
for Neel order and TV = 24. There are two and four sign 
changes as expected. N = An chains have n points 7™ 
with broken Ci or P symmetry. Additional degeneracy, 
if any, is exclusively due the broken Ci and P symmetry 
that indicates VC order. Returning to Figs. [5] and [31 we 
see that lines 7„i — ^^(A) through the mth degeneracy, 
either Ci or P, of finite systems partition parameter 
space into regions in which C{p) changes sign 2m — 2 
and 2m times. 

The gs energy per site, eg (a. A), is related to spin cor- 
relation functions of first and second neighbors 

eo(a. A) = -C^{1) - AC,(1) + aCi(2) + aAC,(2)(10) 

Degenerate eo(ck,A) at 7„i = a„i(A) is achieved with 
unequal C(l; 7™), C(2; 7,„) in Figs. [7]and[5] Finite-size 
effects generate small discontinuities in S{a, A) at all 7™ 
with unequal C{p; 7m). It is straightforward to count the 
number of Ci or P degeneracies in finite systems. We 
cannot predict the number of VC degeneracy, the stars 
in Fig. [2]and|3l and a finer grid may reveal additional gs 
degeneracy leading to VC order. 

In contrast to P or d degeneracy, degenerate eg (a. A) 
under reversal of both Ci and P lead to equal spin corre- 
lation functions as shown in Fig. ^ upper panel, for N 

20, A = 0.75, a = -0.338 and -0.357. The first point 
has four sign changes, in accord with two P degeneracies 
at smaller —a. There is a Ci degeneracy between the 
two VC points, which accounts for six sign changes at 
—a — 0.357. We always find equal C(p;7,„) at VC de- 
generacies within our 3-4 digit numerical accuracy, except 
for = 4n -I- 2 systems with —a > 0.5 where equality is 
limited to two digits and there are pronounced An, An + 2 
effects. Equal C{p; 7^) for all p implies equal S{'jm) ac- 
cording to Eq. ini The VC order parameters Kz{p;jrn), 
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FIG. 8: (Color online) Spin correlation functions C{p) of a 
24-site chain with spin parity degeneracy at 02 = —0.312, 
A = 0.83. 



Eq. m in the lower panel of Fig. [9] has the same period- 
icity but they vanish at points that are shifted by 7r/2. 
It follows from Eqs. [7] and |8] that when 7™ refers to VC 
degeneracy, 

{S^S^_^_.p) = -C_i(p,7m) -iKz{p,'y,n) (H) 

VC order parameters are closely related to transverse 
spin correlation functions in chains with axial anisotropy. 

A simple classical picture captures the principal 
features of quantum spins with VC order ±«;^(p). 
Right and left-handed helices are doubly degeneracy 
with equal energy eo(7m) and spin S'(7m)- H{a,A) 
has axial symmetry in spin space and spins in the 
xy plane for A < 1. We picture the two helices as 
(a;„,±2/„) = (cosn0, zbsmnc/)) with pitch angle (f> that 
generates a specified number of sign changes between 
n = 1 and N. Inversion gives n ^ —n and interchanges 
the helices. Spin parity reverses the z and y components 
of spin, but not Sx = {s'^ + s~)/2, and also interchanges 
the helices. Equal C{p,jm) follow from projecting the 
helices on a plane that includes the z axis, as illustrated 
by 2{xnXn+p)n = 2(j/„2/„+p)„ = where < .. >„ 

indicates an average over n. The VC order parameters 
Kz{p,Jm) go as ±2{xnyn+p)n = ±sinp(j) , shifted by 7r/2. 

When H{a,A) in Eq. [T]has classical spins, the gs is 
easily shown by energy minimization to be a spiral with 
pitch angle (p that depends only on a, 

cos(p = -J1/AJ2 = -1/4q! (12) 
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FIG. 9: (Color online) Spin correlation functions C(p) and 
order parameters Kz{p), Eq. [71 of a 20-site chain with VC 
degeneracies at A = 0.75 and a = —0.338 and -0.357. 



A finite chain with PBC has N(p = 27rm. The critical 
point — 4q! — 1 corresponds to </? = 0, the ferromagnetic 
gs. Finite N ~ An limits the remaining angles to n—1 val- 
ues < If < tt/2. Spiral spin correlation functions appear 
automatically in DMRG treatments j2|| iSlj of isotropic 
spin-1/2 AF/AF chains. This follows directly from a sin- 
glet gs and Eq. [8]with 2C^{p) = C±{p), ahhough DMRG 
is not accurate for the sum over all spin correlation func- 
tions. A spiral interpretation of C{p,"fm) holds in finite 
F/AF chains at any 7^ or, indeed, at any (a. A). There 
are crucial differences, however, between helices and spi- 
rals. Left and right-handed helices are degenerate at 7™ 
with reversed C; and P symmetry. Spirals are associated 
with Ci — ±1 and BOW phases, but are not degenerate 
in ±ip in Eq. 1121 The pitch angle ip — (j) oi spiral or 
helices of classical spins follows immediately from PBC. 
Table H] lists degenerate 7™ at A = 1 for quantum spins 
that can be converted to angles ipm using Eq. [121 The 
proper number of 7^ is found for each iV, but the actual 
values differ considerably for —a > 0.4. Moreover, the 
points 7m = am(A) vary substantially with anisotropy. 



IV. DISCUSSION 

Quite unusually for model Hamiltonians, finite 
spin-1/2 chains with iJ(a, A) in Eq. [H periodic bound- 
ary conditions and frustrated exchanges Ji < and 
J2 — aJi > have doubly degenerate ground state 
(gs) at many parameter values 7m = (am,Am). There 
are three kinds of gs degeneracy: inversion symmetry 
Ci = ±1 or spin parity P = ±1 or reversal of both Ci 
and P. We have exploited degeneracy for finite N using 



FIG. 10; (Color online) Approximate gs phase diagram of the 
anisotropic F/AF chain H{a,A), Eq. [H The phase bound- 
ary between high and low spin is an excited-state crossover 
discussed in the text. The degeneracy line ai(A) has Neel 
order between and As, dimer order for A < Aa- Dashed 
lines enclose regions with multiple degeneracies (inversion, 
spin parity, both) where phases are not assigned. 



exact diagonalization to construct broken-symmetry 
gs according to Eq. [H Exact gs make it possible 
to compute order parameters for broken symmetry, 
spin correlation functions C{p,jm) in Eq. [51 and the 
normalized expectation value 5'(q!, A) of the total spin 
in Eq. [31 Access to exact gs properties compensates 
to some extent for the inherent limitations of finite-size 
approaches to infinite chains. 

We have focused on gs degeneracy to construct and 
characterized broken-symmetry states. Quantum phase 
diagram derived from finite systems invoke excited-state 
degeneracy [3, ll^]. Our finite-size results are summa- 
rized in Fig. [TUl The indicated phase boundary is based 
on the excited-state degeneracy between first excited 
state in Sz = and the Sz = ±1 gs up to TV — 28, as 
discussed in related ID systems 0, [27l - l29| that conserve 
S. The N dependence of the excited-state degeneracy 
is comparably weak and independent of A. The order 
parameter becomes nonzero and an energy gap opens 
very slowly 0, H^l at a Kosterlitz-Thouless transition. 

The high-spin phase [13] is a ferromagnet or a 
Tomonaga-Luttinger liquid (TLL). We find the tran- 
sition between high and low spin to be first order for 
Ab > 0.95, continuous for A < A^. Degeneracy along 
71 — ai(A) is exact at A = 1 or and depends very 
weakly on N in between; as seen in Fig. \E\ there is Neel 
order along along the line between Aa and A b , dimer or 
BOW order for A < A^. An almost identical boundary 
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is shown in Fig. |4] of ref. 17a or Fig. 2b of ref. 17b, 
based on a = J2/I Ji| as a perturbation to an exact field 
theory at A = 0. As noted in connection with Fig. [SJ 
the order parameters of the dimer and Neel phases also 
agree weh. 

Multiple degeneracy is limited to intermediate 
anisotropy A and frustration ratio a close to ai{A). 
The dashed lines in Fig. [10] enclose the parameter space 
in which we find all three kinds of gs degeneracy. We 
cannot assign phases in this region. The degeneracies 
in Figs. [5] and [3] do not evolve systematically with N. 
Finite-size results merely place restrictions on phase 
boundaries. The isotropic (A — 1) chain is limited to 
a BOW or dimer phase for —a < 1/4, consistent with 
exclusively Ci degeneracy in Table HI 

Broken Ci symmetry leading to a dimer or BOW 
phase dominates on the low-spin side. We find i?z(7m), 
the z part of the order parameter in Eq. [5l to change 
sign with decreasing A, from singlet-type pairing at 
A 1 where the gs is a singlet to triplet-type pairing 
that is exact atA = 0,Q: = — 1/2. A dimer triplet phase 
appears for small A in Fig. 2b of ref. 17b, separated by 
a VC phase from a dimer singlet phase. We do not find 
evidence for VC order outside the dashed lines in Fig. [TUl 
and consider the sign of i?z(7m) to be incidental in the 
dimer phase. The small region of multiple degeneracies 
in Fig. [TU] is the major difference with the extensive VC 
regions in some [17, 19, 20, 22] zero-field phase diagrams, 
but not in othersflil. [23{ . 

Finite-size results complement approximate treatment 
of extended systems. DMRG grows a discrete extended 
chain while field theory deals with a continuum version 
of the chain. DMRG [H [3l| results for F/AF or 
for AF/AF spin- 1/2 chains are for open boundary 
conditions and even N that increases by two or four 
sites per step. Even chains have inversion symmetry at 
the center of the middle bond, but not Ci at any site, 
and Ci is relevant for gs degeneracy. Exact treatment 
of a half-filled band of N free electrons with o pen 
boundary conidtions leads to a nondegenerate gs [28l l3l| 
with a BOW whose amplitude B decreases as 1/N and 
excitation energies of order 1/iV. DMRG works well 
for model parameters leading to substantial B > 0.01 
such as the broad peak around a w —0.4 for isotropic 
F/AF chains in Fig. [H But DMRG fails in extended ID 
systems where, as discussed by Affleck and Lieb [l^, the 
distinction between gs degeneracy and 1/A'^ excitation 
energies has to be considered. 

Turning to field theory, we note that a continuum 
model is an approximation, that additional approxi- 
mations are typically needed, and that different field 



theories are possible for H{a,A). Two field theories 
[2^ . [2^ for isotropic AF/AF chains are not limited to 
Ji > 0. The anisotropic F/AF chain has been treated 
with first-order corrections in J2/J1 starting 17[ with 
exact field theory at J2 = and also as bosonization [l^ 
in the opposite limit of | Ji | ^ J2 ■ The merits of field 
theory are beyond the scope of this paper. A continuum 
models suppresses the important distinction between 
inversion symmetry at sites and at the centers of bonds. 

We have discussed frustrated chains with axial 
anisotropy for parameters leading to gs in the Sz — 
sector. The energy spectrum of isotropic (A = 1) chains 
in a magnetic field h = q^bH that defines the z axis 
is simply the zero-field energy plus the Zeeman energy 
hm, with TO = ±1,±2, ... ± S for states with spin 
S. The system has axial symmetry for h > Q but spin 
parity is no longer conserved. In strong fields, the gs 
has multipolar phases near the boundary between high 
and low spin [10|. The general problem of H{a,A) plus 
a static field is considerably more complicated because 
Sz is conserved only when h is along the unique axis. 
Otherwise, the gs has to be found by separately as a 
fimction of h. 

In summary, finite F/AF models H{a,A), Eq. [U 
with frustration ratio a — J2/J1 and axial anisotropy 
< A < 1 have doubly degenerate gs at multiple points 
7m = {am, Am), all in the Sz — sector. Exact gs make 
it possible to compute order parameters at 7m , spin corre- 
lation functions and the spin per site, ^(q;, A). The tran- 
sition from high to low spin with increasing —a is first 
order for A > 0.95, continuous for A < 0.95, with stag- 
gered magnetization or Neel order for 0.72 < A < 0.95 
and dimer order for stronger anisotropy A < 0.72. When 
both Ci and P symmetry are broken, vector chiral (VC) 
order leads to identical spin correlation functions and 
S{a, A) that is interpreted as right and left-handed he- 
lices along the unique axis with spins in the xy plane. 
Finite-size effects are typically quite small, small enough 
to discuss infinite chains. Exact finite-size results are 
consistent with with many-body treatments aside from 
limiting a possible VC phase to intermediate a and A 
where S{a, A) varies rapidly but continuously. 
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